GAMOW VECTORS IN A PERIODICALLY PERTURBED 
QUANTUM SYSTEM 



M. HUANG 



Abstract. We analyze the behavior of the wave function ip{x,t) for one di- 
mensional time-dependent Hamiltonian H = —d^ ±25{x){l + 2rcosujt) where 
■>p{x, 0) is compactly supported. 

We show that ip(x,t) has a Borel summable expansion containing finitely 
many terms of the form E^^L-oo e'^^V-^fc+""'l^l Afc,„e-^**+"'^", where 
represents the associated resonance. This expression defines Gamow vectors 
and resonances in a rigorous and physically relevant way for all frequencies 
and amplitudes in a time-dependent model. 

For small amplitude (|r| <S 1) there is one resonance for generic initial 
conditions. We calculate the position of the resonance and discuss its physical 
meaning as related to multiphoton ionization. We give qualitative theoretical 
results as well as numerical calculations in the general case. 



1. Introduction 

Gamow vectors and resonances, introduced by Gamow to describe a-decay (cf. 
[T]), are very important mathematical tools in the study of metastable (or qua- 
sistable) states in quantum mechanics (cf. [2])- The decay states described by 
Gamow vectors are also linked to the Fermi-Dirac golden rule (cf. [9]). There are 
numerous definitions of resonances and resonant states, using the scattering matrix, 
rigged Hilbert spaces, Green's function, etc. (cf. [Hllini and the references therein) 
These definitions rely on the time-independent Schrodinger equation, though they 
may be extended to time-dependent settings in a perturbative regime (cf. [31 [T3]). 

In a recent paper [15], the author and his collaborator gave a rigorous defini- 
tion of Gamow vectors and resonances for compactly supported time-independent 
potentials in one dimension, using Borel summation (for a detailed description of 
Borel summation, see [151 [5] )■ In this paper, we study the resonances associated 
to a time-dependent periodic potential. In our case, the Gamow vector is of the 
form of the so-called Floquet ansatz (cf. [3]). Our result holds for all amplitudes 
and frequencies of the time-dependent field. In the case of small amplitude or high 
frequency, we calculate the resonances asymptotically, and the real part of the res- 
onances measures the ionization rate. In this sense, our paper extends the results 
of [51 116j . As we will see, time dependency introduces new subtleties and complex 
phenomena. 

2. Setting and Main Results 
We consider the time-dependent one-dimensional Schrodinger equation 
d h d'^ 
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where the potential V{x, t) is a deha function potential well or barrier with a time- 
periodic perturbation. In this paper, we consider two simple but illuminating cases: 

(1) delta potential well V{x,t) = -2A5{x){l + 2rcoswt) 

(2) delta potential barrier V{x,t) = 2A5{x){l + 2rcosLot) 

Here A > Q represents the strength of the potential, r represents the relative 
amplitude of the perturbation and lo the frequency. Without loss of generality we 
take r > OjCJ > 0. We further assume the initial wave function V'o(2;) := ^l^{x,Q) is 
compactly supported and on its support. 

We first normalize the equation by changing variables x 2m3^' * ~^ ImA^^^ ^ ~^ 
"^^f LO. Note that this is more than using atomic units since we also used the special 
property of the delta function = b(x)IA. The equation becomes 



(1) i—%l){x,t) = -^V(a;,t)T2(5(a;)(l + rcoswt)?/'(a;,t) 

(where "-" corresponds to the delta potential well and "+" corresponds to the 
barrier) We shall focus on the delta potential well and analyze in detail the behavior 
of the wave function as well as the resonances of the system for all amplitudes and 
frequencies. The analysis of the delta potential barrier is very similar and we will 
give the results in Section [4] without detailed proofs. 

Theorem 1. Assume the initial wave function 7/1(2;, 0) is compactly supported and 
on its support, then we have for all t > 



k—1 n— — 00 

1 ^ ne^^GO 1 ne^^ 00 

— V / g»^/V-9+"'^»kl+«c.*t-9t^^(_^)^g / F{x,-q)e-'"dq 



27ri 

n— — 00 



where Afc + mui are resonances of the system (Re(Afe) > 0), ip a ramified analytic 
function with square root branch points at every nuji (n g "L), and F an explicit 
function with y/pF{p) analytic in ^Jp. 9 is a small angle chosen to ensure that no 
resonance lies on the path of integration. 

Moreover, the coefficients Ak^n satisfy the recurrence relation 

(2) (^\/~i\/i + nui - Afc - 1^ Au^n = rAk^n-i + rAk,n+i 

and ip{x,t) has the Borel summable representation 



V'(x,t) = i-"^r 



' n=-oofc=0 
K 00 



,-3/2 

e 

fc— 1 n— — 00 

Corollary 1. For 1 ^ fc ^ K , the Gamow vector term 



-\kt+nuit 
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is a generalized eigenvector of the Hamiltonian, in the sense that it solves but 
grows exponentially (in a prescribed fashion) for large \x\. 

Proposition 2. For small r there is only one array of resonances, i.e. K ^ 1. 
The asymptotic position of the array of resonances and a similar result for large lo 
are given in Section\3.4\ 



In the above formulas the branch of the square root is chosen to be the usual 
one: arg(z) e (— 7r,7r] and a,Tg{^/z) € (— fif]- We refer to this choice of branch 
when we use the phrase "usual (choice of) branch" in this paper. 

For small r we calculate asymptotically the position of the resonance, which is 
related to the ionization rate. For generic r we will give numerical results showing 
that the Gamow vector terms exist for some but not all r, and we plot the graph 
of the positions of resonances with different amplitudes (see Section 14]). 

Remark 1 . Theorem [I] and its corollaries generalize to the case where 
V{x,t) = T'^A6{x) 1 + 2^ (rfe cos kujt + Sfc sin kujt) 

3. Proof of Main Results 

3.1. Integral reformulation of the equation. We first consider the Laplace 
transform in t 

/•oo 

■4!{x,p) = / e~P*'ip{x,t)dt 







The existence of this Laplace transform (for Re(p) > 0) follows from the existence 
of a strongly differentiable unitary propagator (see Theorem X.71, [6] v. 2 pp 290, 
see also [7], [H] and [H]). As we will see, Theorem 1 follows from analyzing the 
singularities (poles and branch points) of the analytic continuation of ^j{x,p). 

Performing this Laplace transform on ([1]), we obtain 



(3) ipipix^p) - iipo{x) = 

— -—^ip{x,p) — 25{x)'4>{x,p) — 2r6{x)ip{x,p — iui) — 2r5{x)ip{x,p + iuj) 
ox-^ 

We then rewrite the above ordinary differential equation as an integral equation 

by inverting the operator + ip. We have 

ox^ 



^(x,p) = — / e' •^''g{s)ds — / e' '^'g{s)ds 

where 

g{x) — iijja{x) — 2S{x)'ip{x,p) — 2r5(x)ilj{x,p — iuj) — 2rS{x)ip{x,p + ito) 

Recalling that 5{x)f{x)dx = /(O), we simplify the above integral equation 
and obtain 
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-0(0, p) + r'(/;(0,p — iuj) + rip{0,p + iuj) 



VP 

Letting x = we get an equation for ip(0,p) 



(5) V(0,P) - / ^yMs)ds -TT^ e-^' ^^Ms)ds 



+ — ^-0(0, p) + r0(O,p — iw) + r?/;(0,p + iuj)^ 



VP 

which impHes 

(6) (0(O,p) + r0(O,p - icj) + r'tp{0,p + iuj)J = 

Substituting ([6|) in ([4]) we get 

(7) 0(x,p) = e*'''v15kl^(o,p) + _ e''''v^l-l/(0,p) 
where 

^Vp J+oo ^Vp j-oo 

Equation ([7|) indicates that the analytic continuation of as weh as its 

singularities, follows naturally from that of 0(0, p), so it suffices to analyze 0(0, p) 
using the recurrence relation ([5]). Later we will perform the inverse Laplace trans- 
form on ■0(x,p), justified by estimating 0(0, p) and f{x,p) for large p. We will 
then deform the contour of the Bromwich integral, which yields the expression in 
Theorem [TJ It is worth noting that to deform the contour it suffices to place a 
branch cut of the square root in the left half complex plane, while to analyze the 
singularities of 0(0, p) we need to consider a larger region in the Riemann surface. 
Some delicate points of the analysis stems from the complexity of the Riemann 
surface, since, as we will see, ip{0,p) has infinitely many branch points and there 
appears to be a barrier of singularities on the non-principal Riemann sheet. 

3.2. Recurrence relation and analyticity of 0. We rewrite the recurrence re- 
lation ([5]) as 

(V^Vp^ ■0(0, p) = r0(O,p - iuj) + r0(O,p + iuj) + V~iVpf{^-,P) 

We will show that f{0,p) — ^ ^^^^ -\- Q ( as p — > cx) in any direction in the 

p yp'^'^y 

right half complex plane (see Section [3.6|1 . It is not a priori clear that 0(0, p) has 
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an inverse Laplace transform. We thus let ipip) — 1^(0, p) — f{0,p). The recurrence 
relation for ijj is 

(8) (y^s/p - l) i>{p) = r^pip - iw) + r'4){p + ito) + (1 + 2r)/(0,p) 

It is convenient to write the recurrence relation in a difference equation form. 
Denotingp = i+inu>+z , yn{z) = ^{i+intu+z), and fn{z) = {l+2r)f{0,i+inLU+z), 
we have 

(9) (y^y/ i + inuj + z- 1^ yn{z) = ryn-i{z) + ryn+i{z) + fn{z) 
The associated homogeneous equation is of course 

(10) (y^V i + inuj + z - ij yn{z) = rt/„_i(z) + ryn+i{z) 

Let Zq be a branch point closest to 0, that is, a point on the imaginary axis 
satisfying —z^i = inf„{|l+nti;|} (note that l^ol ^ and let no be the correspond- 
ing n. Since clearly y-niz) — y„+i(z — iuj) = ?y„_i(z + iuj), it suffices to consider 
Im(z) e (—^uj, ^uj) for the usual branch. In general, if we make a branch cut 
at {e^^oo,zo) {cos6 ^ 0) we consider the strip-shaped region {|Im(z) — psinffj < 
|w,Rc(z) =p cos 0,p e R}. 

To analytically continue y := {yn}, we consider the Hilbert space Ti defined by 



\n\'/')\xJ' 



n= — oo 

and the operator Cm :Tl ^ H 

s /X {I + mVi)yn{z) + ryn-i{z) + ryn+i{z) , 

{Cmy)n{z) = ^ 7=r (m £ Z+) 

It is easy to see that Cm is entire in r and analytic in ^z — zq in the region 
Re{z)>-m\lm{z)ei-lu;,lio). 

Lemma 3. Cm is a compact operator for any choice of branch. 

Proof. For arbitrarily large € N, we consider the finite rank operator 'Dm,N '■ 

' iCmy)n \n\<N 

otherwise 



n 

It is easy to check that 



\\Cm-Vm,N\\=0{.N-'''') 

Therefore Cm, being the limit of finite rank operators in operator norm, is com- 
pact. 

□ 

Lemma 4. The equation 

{y^Wi + inu) + z- 1^ yn{z) = ry„-i{z) + ryn+i{z) + gn{z) 
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has a unique solution in Ti for |Re(2)| > (2r + 1)^, for all g G ?i. In particular, ([9|) 
has a unique solution and ()10|) has only the trivial solution y — 0. The conclusion 
holds as well if z ^ Q and r is sufficiently small. Furthermore, for large |Re(2;)| we 
have |y| = 0{\Re{z)\~'^/'^\^\) where |x| sup„ |x„|. 

Proof. Note that under the assumptions above, the norm of the hnear operator 

s -.n-^n 

\ ( \ - ryr,^i{z) + ry„+i{z) 



(^^/~t\/ i + inuj + z — l) 



is smaller than 1, since | \f^i\Ji + inuj + z — l| \^i + ini^ + z\ — 1 ^ ^|Re(z)| 
1 > 2r. We then have 



y = 



(^/— 1\/ i + inuj + z — l) 

□ 



Proposition 5. For every r € C, there are at most finitely many z = z\,...,zi^^ 
for which the homogeneous equation (jlOp has a nonzero solution y inTL. For all 
other z, there exists a unique solution to (jO]). The function \/ z — z^y is analytic in 
both '^J z — zq and r, and it can be analytically continued on the Riemann surface of 
\/i + incu + z to argz G (— 37r/2, 37r/2). (in other words, one can rotate the branch 
cut in the left half complex plane) Moreover, zi, zi^ are either poles ( in \/ z — zq) 
or removable singularities of y, and yn {n 7^ no) is analytic in -y/z — zo when z is 
close to Zq. 



Proof. We consider the equation 



1 



(a/^V* + inu) + z + m\/I) 



Since Cm is compact, analytic in both r and ^/z~^^~zo, and invertible for |Re(2)| > 
(2r + 1)2, it follows from the analytic Fredholm alternative (see [B| Vol 1, Theorem 
VI.14, pp. 201) that the proposition is true for every y^™' (note that the solu- 
tion of the inhomogeneous equation exists for |Re(z)| > (2r + 1)^, thus there can 
only be finitely many isolated singularities). Uniqueness of the solution implies 
yM = y["i+i] for all r G C,Re(z) > — m^. Thus we naturally define the analytic 
continuation of the solution to be y := y[™l . Analytic continuation on the Riemann 
surface follows from the fact that for fixed r, z {z not on the branch cut) slightly 
rotating the branch cut does not change the value of + inuj + z for any n e Z. 
Uniqueness of the solution thus ensures y also remains unchanged. 

Assume yn{z) ~ 6„(z — zq)^^^^ as z ^ zq. It is easy to see from © that 

(y^Vi + inuj + z - ij bn = rbn-i + rbn+i in ^ no) 

W~^i\Ji + mow + zo - l) = rb^o-i + rb^o+i - (l/2 + r)i^/^ / il;o{x)dx 
The unique solution of this recurrence relation is obviously 

/oo 
V'o(x)dx,6„ = 0(n^no) 
-CXD 

□ 
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Corollary 6. For every r E C, (0) has a unique solution tp. y/pip meromorphic 
in p with square root branches at every inoj (n £ Ti) and poles at {pk + inuj} 
(k ^ l,2.../r,n e Z). 

Proof. In order to recover p — i inuj + z from the solution to ([5]), we only need 
to show yn{z) = 2/n=Fi(-z ± 1^*)- To this end, note that by © we have 



(11) {^\/~-i\/ i + inoj + z — \ j yn^i{z ± Loi) 

= ryn^i^i{z ± uji) + ryn^i+i{z ±uji) + /„=pi(^ ± uji) 

which is the same equation as ^ since fn^iiz ± uji) = /„. 

Thus, uniqueness of the solution (Proposition [5]) implies yn{z) = yn^iiziuji). 
Note that we need to choose the same branch for all \/i + inoj + z. □ 

We conclude this section with a few observations about the positions of the poles 
of V', including the well-known result of complete ionization (see [5l \7\ [16]). 

Proposition 7. For r > 0, y has no pole on the imaginary axis or the right half 
complex plane, with the usual choice of branch. 

Proof. In view of Proposition O we only need to show the homogeneous equation 
(jlOP has no nonzero solution in Ti.. Multiplying (jlOp by y^{z) and summing in n 
we get 



oo 

E 

n— — oo n— — oo 



(y^Vi + inuj + z-l^ \y„\'^ = 2r ^ Re{yn-iy„) 



which implies 



\/—i\/ i + inu) + z \ y„ 



must be real. 

If Re(2:) ^ then lm{^/^\/i + inuj + z) ^ for all n and Im(\/^V i + inoj + z) < 
for ah n < -(1 + \z\)/uj. Thus yn = for all n < -(1 + \z\)/uj and ^ implies 

y = o. □ 

Proposition 8. For r > 0, y has no pole on the imaginary axis for any choice of 
branch. 

Proof. Similar to the above. Note that Re(z) — implies lm{^/^\/i + inuj + z) = 
for all n > —(1 + lm{z))/uj and Im(V^V* + inoj + z) has the same sign (and 
nonzero) for all n < —(1 + Im(z))/w. □ 

Proposition 9. Solutions of the homogeneous equation (llOp exist in negative con- 
jugate pairs, in the sense that if zi is a pole of ip, then — zT is also a pole (with a 
different choice of branch, see proof and comments below). 



Proof. Simply note that {—iy/^\/i + inuj + z = {—i^/'^y/i + inuj — z if we choose 
the branches in such a way that in the upper half complex plane the two square 
roots are the same, while in the lower half plane they are opposite. □ 

In view of the above propositions, we will concentrate our study of resonances 
on the left half complex plane. The author believes that the imaginary line on 
the non-principal Riemann surface is a singularity barrier, and the Proposition [5] 
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provides a pseudo-analytic continuation across the barrier. We will not discuss the 
details in this paper. 

3.3. The homogeneous equation. As we mentioned in the introduction, poles of 
y in the left half complex plane correspond to resonances of the system. According 
to Proposition m finding these poles is essentially the same as finding solutions to 
the homogeneous equation (IIOII in H. 

Lemma 10. Assume the nonzero vector u = {un} satisfies the homogeneous re- 
currence relation ()10p . and that 

oo 

^(l + |n|3/2)|z,„|2 <oo 
Assume also that the nonzero vector v — satisfies (fTO]) and 

J2 il + \n\'/')K\'<oo 

n— — oo 

Then the homogeneous equation ()10|) has a nonzero solution in Ti. if and only if 
the discrete Wronskian W := u„Wn+i — w„u„+i — 0. The solution, if it exists, is a 
constant multiple of u (or equivalently vj. 

Proof. If r = the lemma is trivial. Assume r > 0. We first note that the recurrence 
relation ([T0|) implies 

(1) is independent of n. 

(2) for any n and any nonzero vector x satisfying that recurrence relation, we 
have |x„|2 + \xn+i\^ ^ 0, \xn\^ + |a;„+2p 7^ (n 7^ -1). 

Now assume W = 0. Since v is nonzero, there exists m for which v,„ 7^ 0. Thus 
we have Um±i ~ {um/vm)vm±i- Since u 7^ we must have Um 7^ 0, for otherwise 
Um±i = Um = 0. If Vm±i = then Um±i = 0, which implies \um - (w,„/um)w„ip + 
\um±i - ium/vm)vm±i\'^ = 0, meaning u = (um/um)v. If Vm±i 7^ then Um±l ^ 0, 
which inductively implies again u — (um/t'm)v. Therefore u solves ^TO\i in H. 

If W then clearly u and v are the two linearly independent solutions of the 
second order difference equation (jlOp . Furthermore, we have liminf„<o |u„| > 
and liminf„>o |wn| > 0, since limsup„>Q |u„| < co7isi.|n|^'^/^ and limsup„<o \^n\ < 
const.\n\~^^^ but UnVn+i — VnUn+i is a nonzero constant. Therefore no nonzero 
linear combination of u and v can be in 7i. Since a second order difference equation 
cannot have any other solution, there is no nonzero solution of (jlOp in Ti. □ 

We now give a constructive description of u and v. For convenience let hn{z) = 
[^/—i^/ i + inuo + z — l). We choose ni,2 € Z so that |/i„| > 2\r\ for all n ^ ni > 
and n ^ 712 < 0. Let I be the identity operator. We define Tiia by 



oo 
n—ni 

712 

n— — oo 
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Proposition 11. There exist u and v, analytic in r and ramified analytic in z, 
satisfying the conditions described in Lemma \7R Moreover, u{z±uji) — const.u(z) 
and v(z ± Loi) — const.v{z). 



Proof Let Ti : Til ^ 



Let a= {r/hn,, 0,0...) £ Hi. 
The equation 



has a unique solution 



r 

-j—{yn-i + Vn+i) n> ni 



r 

yn+i n = ni 



hn 



Tiu 



u = (X - Tiy^a = a + Tia + T^a... 
since clearly ||7i|| < 1. It is easy to see that u satisfies 

r 

— (u„_i+u„+i) n>ni 
r 

— (u„+i + l) n = ni 
. hn 

Thus the recurrence relation 

hn 

Un = Un+1 - Un+2 

r 

extends u to a solution of the homogeneous equation ([10]). In particular Um-i = 1- 
This solution u is analytic in r and z (ramified) locally since 71 and /i„ are analytic 
in r and z (ramified), and the uniform limit of analytic functions is analytic. As r 
or |Im(z)| increases we may analytically continue u by considering some > ui so 
that \hn\ > 2|r| for all n ^ 713. Using the same procedure as we did for m we get 
u. It is easy to see that u = Un-^u for they both satisfy the contractive recurrence 
relation (in the sup norm) 

r 

— (u„_l + Un+l) n > 713 

r 

— {un+i+Uns) n = na 
Note that this implies u„ 7^ for large n. 

The analytic continuation of u is, up to a scalar multiple, periodic in z. Note 
that u^{z) = u{z ± oji) satisfies (for large 77,3 > ni) 

-K^-i+"n+i) n>ni 




10 



M. HUANG 



while u satisfies 



Un±l 



hn±l 

r 



Thus u(z ± Loi) 



. hn±l 
U„3 (z ± LOi) 



{Un±l-l + U„±l+l) n > 723 
(Wn±l+1 + Mn3±l) n = 



■"n3±l(^) 



u(z). 



The construction of v is very similar, namely w — {I — T2) where 7^ : TL2 



— {yn-i + Un+i) n <n2 

fin 



n = n2 



and b = (■■■,0, 0, r/hn^)- 



□ 



Proposition 12. 14^ is analytic in r and ramified analytic in z. Moreover, W{z) = 
if and only if W{z ± Loi) — 0. 

Proof. The first part is obvious. The second part follows from the relation u(z ± 
Loi) — "'^ ^ -— ^u(z) (see the proof of the previous proposition) and the fact that 



Un3±l(^) 



□ 



Remark 2. Another way of constructing u and v is by using continued fractions, 
see [S] . The continued fraction expression is slightly simpler in this particular case, 
but our iteration method can be easily generalized to trigonometric polynomial po- 
tentials mentioned in section\^ 

3.4. Resonance for small r. We assume r > and analyze the resonances of the 
system for small r (relative to uj) by locating zeros of W , in view of Lemma [TOl 
Since we will need to consider different branch choices, we write for convenience 
hn{z) = {{—iy^'^\/i + inuj + z — 1) where the power 1/2 always indicates the usual 
choice of branch. 

Lemma 13. For every choice of branch, there exists a constant c so that when 
uj > c{r + r'^), we have |/i„| > 2r for all n ^ 0. 

Proof. Recall that for a branch cut at (e*^oo,zo) (cos^ 7^ 0), we consider the strip- 
shaped region := {|Im(z) — psin6'| < |u;,Re(z) —pcoa9,p£R}. It is easy to see 

that ci := inf„-^o zefit I inl > 0. Therefore |/i„(2;)| = |(— + inuj + z—l\ = 

' LU 

linuj + z\ linu) + z\ cioj „ ^ „ _ 

:^ ' ' :^ > 2r if > 2r + 2^?. 



|V« + into + z + Vi\ ^\inu! + z\ + 2 ^/cTuj + ' 



Note that 



is an increasing function for a: > 0. 



□ 



Proposition 14. For small r, there is a unique nonzero solution of the homoge- 
neous equation U0\) in the left half complex plane with the usual choice of branch. 
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Moreover, the solution satisfies 

( 2i 2i 



y(l+C^)l/2_l i-l/2^(l_^)i, 

where a{r) is analytic in r and (t(0) = 0. 



cr(r) r 



2 



Proof. Wc choose ni ~ l,n2 = —1 to construct u and v. Thus uq — — \ and 
W = v\ — u\. We calculate by iterations 

3 5 

^'1 = W_l = -2-^^ -g i?2 

r /ll /i^/i2 h'^ih-2 h\^^ h\^^ 

where i?i,2 are bounded for w > c(r+r^). Note that |/io(2)| = Wi + z—i^/'^\ ^ |z|/2 
r 



and 



/ir: 



^ — r for all n ^Q. 

ClUJ 



Now, if cj is fixed and r is small, W = implies ho{z) — 0{r^). Hence we must 
have z = O(r^). In addition, we need to make the choice of branch so that ^/i is in 
the first quadrant. Thus we let z = (oq + (T)r^ where a = o(l), and we see that 

oil)) 




' (1 + Cj)l/2 _ 1 i-l/2^(l_^)i_ 1 

For small r, W is clearly analytic in both r and a. Since the value of W depends 
only on ■ \z — inul < 2aor^}, there are exactly two different W with different 

choices of branch, namely Wi : Re(\/i) > 0, Re(\/— i) > and W2 : Re(\/I) > 
0, Ke{\/^) < 0. However, according to Proposition [5] and Proposition [51 they 
are in fact negative conjugates to each other, and only one will be in the left half 
complex plane. We thus take W — W\ for its branch is consistent with the usual 
branch. 



W 

|r=0.(T=0 — 

r 

d fW\ . i 



It is easy to verify that 



da\r J 2 ^ ^ 

Therefore it follows from the implicit function theorem that the position of the 
zero of W is given by 

( 2i 2i , \ 2 

- - ' + cr(r) 



\^(l+C^)l/2_l i-l/2^(l_^)i 

where (T(r) is analytic in r and (t(0) = 0. 
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a{r) can be found asymptotically by iterating a{r) — as in the standard 
proof of the implicit function theorem. 

Since the usual choice of branch is consistent with W, the zero of W is visible. □ 

Corollary 15. For r small and a; > 1, the position of the resonance satisfies 

- 1 „ 2y/u; + 1 „. 
Ai ^ r^i. 

UJ UJ 

Proof. The corollary follows from the expression of qq with the usual choice of 
branch. The fact that it is indeed a resonance, i.e. a pole of y, will be established 
in the next subsection. □ 

Remark 3. In the case uj ^ 1+r^, an analogous analysis shows that the position 
of the resonance is given by Ai ^ ^ ^7^'' 

Proposition 16. For small r the poles (in one vertical array) of ^ are simple and 
the residues are nonzero for generic f . 

Proof. We note that the order of the pole of (X — Cm)~^ equals the order of the 
corresponding zero of X — Cm, which is a constant by the argument principle (see 
Lemma [18] below), since I — Cm is analytic in z. It is easy to verify that when r = 
the zero of X — Cm is of order one. Thus the poles are simple. 

Let z = G{r) be the continuous functions satisfying W{G{r),r) — 0, G(0) = 0. 
We consider the residue 

P{r) = ^.i yoiC,r)dC 

J\C-G(r)\=e 

Obviously P{r) is analytic in r. For generic f , P(0) 7^ (in which case y can be 
found explicitly). Thus P{r) 7^ for small r. □ 

3.5. Resonances in general. Having analyzed the zeros of W for small r, we 
proceed to consider the case for general r, as well as the poles of y. 

For convenience we study the region r2g ^ := {z : Im(z) G [psind + — ^uj + 
e,psin6l + f + ^uj + e),Re(z) ^pcose,p £ K}n{z : |Re(z)| < (2|r| + 2)^}, the 
branch cut being placed at (e*^ 00,^0) {cos 9 ^ 0). It is easy to see that there is 
exactly one zero and one branch point inside this region for small r (cf. Section 
I3.4p . We note that as long as z is not located on a branch cut, we may rotate the 
cut slightly without changing W . 

Lemma 17. For every r, W has finitely many zeros in 1J| cos9|>C(,>o ^s.e where Cb 
is arbitrary. 

Proof. By Lemma IH there is no zero for |Re(2;)| > (2|r| + 1)^ and the zeros are 
isolated. Since the Riemann surface of the square root has only two sheets and the 
region Qg^f is bounded, W can only have finitely many zeros. □ 

Lemma 18. Assume for some rp and arbitrarily small e > 0, with the branch 
choice arg(z) G {6 — e,9 + 2tt + e) (~2tt < 6 ^ 27r,cos0 0), W has finitely many 
zeros in fie.e- Then the number of zeros remains a constant if r is close to rg. 
Furthermore, each zero moves continuously with respect to r. 

Proof. The lemma follows from standard complex analysis arguments. Suppose 
W{z,ro) has zeros zi^Z2^...Zm inside flo and Zm+i, ...Zm+i on dfto. Since z g flg^e — 
ilefi iff z — iw e ^lefi — ^le,e, we let Zm+k = Zm+k + iuj {1 ^ k ^ I). We may 



GAMOW VECTORS 



13 



choose smaU e > so that W{z,ro) has zeros 2i,Z2,...2m+i hi flo,f_ for arg(z) G 
{6 — e,9 + 2TT + e), and no other zero in ne,2e for arg(z) G {6 — 2e,9 + 2tt + 2e). Let 
< (5 < e be small so that there is at most one zero or branch point inside any 
circle of radius 25, and W{z,ro) is analytic (with a suitable choice of branch) in 
\z — z„\ < 26. Since W is analytic in both z and r, it follows from the argument 
principle that for r very close to tq 

1 r i^W(C,r) 
^T^^ J\C~z^\=s W{C,r) 
Now we consider the compact region f2' := {z : arg(z) e [0 — e, 6'+27r + e]} p| ^e,e\ 
[Jn=i{^ ■ \z - ^"1 < Clearly \W{z,ro)\ > for aU z e Since W is jointly 
uniformly continuous in z and r, we have r)\ > for all z £ fi', r close to rg. 

Thus the number of zeros is locally a constant and they move continuously with 
respect to r. □ 

Proposition 19. For every r there are finitely many zeros of W in any strip 
{z : Im(z) e [z, z + uj), Re(z) G M} for all choices of branch within \ cos9\ > Cb > 0, 
and the position of each zero changes continuously with respect to r. 

Proof. The conclusion follows from Proposition [121 Lemma [17] and [TH Note that 
we may choose 9 arbitrarily thus covering the whole Riemann surface (except for 
the imaginary lines). □ 

As we have shown in Proposition [5] and Lemma 1101 poles of y are located 
where W — Q. We summarize the results as 

Proposition 20. For generic r and f, y(z,r) has finitely many arrays of poles 
for any choice of branch with |cos0| > Cf, > 0. Their residues Ak^n satisfy the 
recurrence relation 

((-O^^^V* + '^'^^ - Ai - 1^ Ak,n = rAk,n-l +rAk,n+l 

and Afc G Ti. 

Proof. The first part is simply a rephrasing of previous results (cf. Proposition 1 19[). 
The recurrence relation for residues follows from the fact that 

Ak,n ^ 77— f yk.n{C,r)dC 

satisfies the homogeneous equation (jlOp since y satisfies ([9|) and 

/„(C,r)dC = 

'|C-G(r)|=e 

The above expression for P„ also implies G H since, by Holder's inequality 

oo oo „ 

J2 (i + H'/')l^fe,nl'^ E + |yn(C,r)|^d|ci 

n— ~oo n— oo ^|C-G{r)|-e 

= (f E (l + H'/')ly"(C,r)pd|CK sup ||y(C,r)||2 <oo 
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the last inequality following from the continuity of y (see also Section 13.61 below) . 

□ 

3.6. Proof of Theorem 1. As we have mentioned before, we will take the inverse 
Laplace transform of ip and deform the contour, collecting contributions from the 
poles in the process. We first provide the necessary estimates. 

Lemma 21. Assume supp-i/'o G [—M,M], then y/pf(x,p), where f{x,p) is as de- 
fined in Section\^ is analytic in ^ with a square root branch at zero. Moreover, 

p 

for large \p\ . 

Proof. By integration by parts we have 

= 5 / ^'^'a{s)ds 

^— v^M^)ds 

1po(x) .3/2 ,,,, , , i"3/2 rO , , 



The lemma then follows. □ 

Lemma 22. '(/'(p) satisfies 

(1) For any compact region i^i G C which does not contain any pole of 'ip{p), we 
have 

00 

sup ^ (1 + |ri|^/^)|-0(p + n[ji)p < 00 

Pe^i n=-oo 

/n particular, 

00 

sup >^ + nwi)! < 00 



- n— — 00 



For ant/ c ^ 0, IV-Cp)! dp < 

(^5^ For |Re(p)| > (2r + 1)^ we have 



4,{p) = p-i/^o (/(o,p)) = o + o (p 



-2 Mi"/^^ 



Note that the p ^1'^ behavior of tp{p) near the origin does not affect the nature 
of these estimates, so we omit further discussions of that special case. 
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Proof. (1) Recall that tp{i + nivi + z) — yn{z) and that y i.e. 

l|yir= E ii + \n\'/')\yn\'<^ 

n— — oo 

Since y is continuous in z on the Riemann surface of the square root, so is ||y||. 
Compactness of fJi then implies sup^gQ^ | |y| | < oo, from which the first part follows. 
The second part follows from the Cauchy-Schwarz inequality 



V sup \ljj{p' 



y (1 + |n|3/2)-l/2(l + |n|3/2)l/2 _^ , 



OO 

^ Ed 

n— — oo 



|3/2w 



1 E + sup + < 



(2) Note that by Fubini's theorem and Cauchy-Schwarz inequality (cf. part (1)) 
we have 



/ ^{p) dp= 



^{c + riLji + si 

„i oo 



ds 



c + nuji + s?) ds 



n— — oo 



^ sup E] + < OO 

pe[c-si,c+si] n=-oo 



(3) The conclusion follows from Lemma 0] and Lemma [^T] 



□ 



Proposition 23. ii{x,t) = 2^ /^^ e*'^^^\/Pl^'l+P*'0(p)dp+ 2^ /^^ eP*f{x,p)dp, where 
the contours Ci^2 0*^6 as shown in Figure\^ and\^ In the process of deforming the 
first contour, we collect contributions from the poles and we slightly rotate the branch 
cut by a small angle 6 if a pole sits on the usual branch cut. 



Proof. We first note that 



sup 

lm(p)^0 



and 



sup 

Im(p)<0,s6l 



J^/'^ ^p+is\x\ 



^ sup 



„ci|a:|^-Im(p) 



where Ci = sup^g^Re + iv)^ < 00. 

Now, by the Bromwich integral formula 



2m 



c+ioo 



e^*ip{x,p)dp 



1 



C+ZOO 



e'"^"^\^\+P'Mp)dp+ — 



c+ioo 



eP'f{x,p)dp 
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By Lemma [2l1 we have 

^T^l J c-iao 



' c—ioo \ P 



2m 



-dp . ^ . 
p Zm 



dp 



Mx) + ^ f - ^) dp=^ i eP'f{x,p)dp 



2m 



2m 



IC2 \ f / JC-2 

As for the first contour, we only need to show that (along both sides of the 
branch cuts) 



E 



n— — 00 



^i^^^ \/ s+nuji\x\ + st+nujit 



ip{s + nuji)ds < 00 







00 ^-ge''' + (n+l)ui 



and if the resonance is visible with the usual (or slightly rotated) branch cut, then 



\Ak,n\ < 00 



n— — 00 

The first two estimates follow from Lemma [221 since 



J2 I e*'^'^''^+^l^l+''*+"'^'*V;(s + nwi)ds 



n— 00 



< sup ^ \'4)(p + nu}i)\ 



-qe 



^ci\x\^\+st^^ < 00 



and 



00 /•-ge''' + (n+l)i^i 



.i'''^^\x\+pt 



i){p)dp 



„ Jo 



ds 



n——oo 



^^c,\x\^\-qcoset I J2 \i;{p + nuji)\] < 00 

The estimates for the resonances follows from proposition [20] and the Cauchy- 
Schwarz inequality. Since G 7i, we have 

Y: \A,,n\^ J2 (l + |r.r)-^/^(l + |nn^/^|A.,n| 



< E (1 

n— — 00 



|3/2n 



-1 £ (l + |n|3/2)|^,,„|2<oo 



n— — 00 
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□ 



Figure 1. Contour Ci 



Figure 2. Contour C2 



Corollary 24. For t > 0, we may further deform the contour Ci to C3 by pushing 
the vertical lines left to infinity. 



Proof. Note that, in the proof of the previous proposition, '''^ gCikl v|s|+st^g 
bounded in Re{q) > and e'^^^'-'^VM-Qcoset ^ g as Re(g) ^ 00. 

Thus we conclude the proof of Theorem [1] by taking the differences between the 
upper and lower branches to deform the contour integrals into line integrals. To 
be exact, if we denote Fs{x, y/p) = f{x,p), (finiVP ^ nuji) = "0(0, p — ncui), then we 
take F{x,p) = Fs(x,y^) - Fs{x,~y^) and ip{p) = ipn{x,y^) - -0„(a;, ^). 
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The last part the theorem follows immediately from Watson's Lemma, since F 
and If are clearly analytic in y/p and has sub-exponential growth as Im(p) —oo 
(see Lemma HH and mi). Note also that ■i^b) (1 + 2r)/(0,p) as p 0. 

Corollary [1] follows from a direct calculation using ([T]) and ([2]). □ 



Figure 3. Contour C3 



4. Further Discussion and Numerical Results 

In this section we study the physical meaning of the resonances, calculate the 
positions of the resonances numerically, and discuss the delta potential barrier. 



4.1. Metastable states and multiphoton ionization. When a resonance is 
close to but not on the imaginary axis, it corresponds to a metastable state of the 
wave function (see [15|). If is not too large, for a moderately long time the wave 
function is governed by the Gamow vector terms whose resonances are closest to 
the imaginary axis . Thus, for a fixed initial wave function, the real part of these 
resonances approximately measure the rate of ionization, that is, the integral of 
l^/jp over a fixed spacial interval as a function of t. 

It has been observed (see 5 ) that the rate of ionization changes rapidly when 
w is approximately equal to an integer fraction of the bound state energy (in our 
case, uj = 1/m, m e N). This phenomenon is related to multiphoton ionization 
(see [5l [11] [T2I [13] and the references therein), a process in which an electron 
escapes from the nucleus by absorbing multiple photons at the same time. Since, 
as we mentioned in the last paragraph, the ionization rate can be measured by the 
position of resonances, we expect a rapid change in the real part of the resonance 
Ai when uj is near 1/m and r is small. 

Proposition 25. For ^^-^ < u; ^ the real part of the resonance is of order 
^2m+2 small r. 
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Figure 4. Real part of the resonance as a function of uj 



Sketch of Proof. Recalling Proposition 1141 wc have z 
It is easy to see that Im(/i„) = O(r^) for n ^ ~1/lo. 



2i 



(l+w)l/2_l ^ 



It can be shown by induction that (7^ v)i is a function of h^2---h_ 



J)T72 ■ 

and 



of order r^+^ . Moreover, {Ti^+^w)i = and {T^^^)i ^ 0. 



Therefore, with the notation z 



a] , we have 



_(l + w)l/2_l I l-(l-w)l/2 

lin(W) = -§Rc(cr)(l + o(l)) - c^r2™+i(l + 0{r)). Thus we must have Re(CT) - 
const. r^"^. 

The above proposition implies that there is indeed a rapid change in the real 
part of the resonance. Here we confirm this result with numerical calculations (see 
Figure m below) and omit further details of the proof. 

4.2. Position of resonance: numerical results. As we have shown in Section 
13. 4[ for small r there is only one resonance in the left half complex plane, for all 
choices of branch. This is, however, not always the case for general r. 

We demonstrate the position of resonances in the left half plane by numerically 
calculating zeros of W for different r. In the graph below we show zeros of W 
plotted with different r and choices of branch, with lu = 2. 

Based on these numerical results, we make the following observations: 

(1) For some values of r, such as those between 0.69 and 1.31, there is no 
visible resonance with the usual choice of branch. In other words, the 
Gamow vector term in Theorem 1 is absent. 

(2) New resonances ("+" marks) emerge as r becomes larger. They can only 
be "born" from the imaginary axis, according to Proposition [T9l 

(3) With any given r, there does not seem to be more than one resonance visible 
with the usual choice of branch. 

(4) Resonances always move upward with increasing r. 

(5) New resonances move farther away from the imaginary axis compared to 
older ones. 

(6) "Old" resonances ("x" marks) do not move arbitrarily close to the imagi- 
nary axis with increasing r. 



4.3. Delta potential barrier. Finally, we briefly discuss the case for the delta 
potential barrier. The corresponding recurrence relation (see ([9])) is 
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r>1.43 



Im(p)' 



0.77<r<1.32 



Retp) 




-1.0 -O.B -C-.4 -0.2 

Figure 5. Position of resonances for different r. Dots are reso- 
nances for the usual branch, and "x" and "+" are those resonances 
continuing on the Riemann surface (they are not visible with the 
usual branch cut). The "x" and "+" curves in the middle are on 
different Riemann sheets. 




(y-i\/p + 1^ '*/'(0,_p) = ril){0,p- iuj) + ril!{0,p + iuj) + \/~~i^f{Q,p) 

With a change of branch ^ ^^"^ changes of variables r — r, / — /, 

the above equation becomes 

[^/~~i^/p- ipiO,p) = ril){0,p- iu) + rV)(0,p + iu) + ^/~^^f{Q,p) 

which is identical to 

Therefore essentially all the theoretical results hold for this case as well. Note, 
however, that for small r there is no resonance with the usual choice of branch 
(which corresponds to a different choice of branch in the potential barrier case, see 
Proposition . 

For larger r, we expect the behavior of the wave function to be qualitatively 
similar to that with a delta potential well, since the contribution from the time- 
independent part will be relatively insignificant compared to the time-dependent 
part. This is confirmed with the graph below plotted for different r and ui = 2. We 
choose the usual branch for simplicity. 

Acknowledgments. The author is grateful to O. Costin who introduced him to 
questions of the type addressed in the paper and gave him many valuable sugges- 
tions. 
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